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The interplay of the electron-phonon interaction and of the double and su- 
per exchange magnetic effects is analyzed in the Lai-xC a^MnO^ perovskites 
with < X < 0.5. By using an analytical variational scheme that allows to 
treat the electron-phonon interaction in a fully quantum manner it is shown 
that in the intermediate electron-phonon coupling regime a charge density in- 
stability occurs near the metal-insulator transition, induced by either temper- 
ature or hole doping. The system segregates in antiferromagnetic or paramag- 
netic small polarons and ferromagnetic large polarons domains characterized 
by different values of the lattice distortions that are satisfactory compared 
with those experimentally observed. 
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Introduction. In recent years the perovskite oxides Lai^^C a^MnO^ liave become the 
focus of the scientific interest owing to the discovery of the colossal magneto resistance effect. 
These materials have been first studied in the 1950's for their peculiar phase diagram where 
is present a strong correlation between magnetization and resistivity.0 In particular, in the 
phase of colossal magneto resistance,il characterized by 0.2 < x < 0.5, these compounds 
exhibit a transition from a paramagnetic insulator to a ferromagnetic metal upon cooling 
that has been generally explained by the double exchange mechanism originally introduced 
by Zener.i The discovery of the colossal magneto resistance phenomenon and the observation, 
by means of the resistivity studies,^ that the tight binding band description of the conduction 
electrons!'! is inadequate and that strong interactions between the electrons and the lattice 
distortions play a non negligible role have aroused a renewed interest for these compounds. 
The relevance of the Jahn- Teller polaron formation, first suggested by Millis,! has been 
shown by giant isotope shift of the Curie temperature,! by measurements of the lattice 
distortions by EXAFS0 and by frequency shifts of the internal phonon modesi and now it 
is generally recognized.! 

Recently, a large amount of experimental data has been collected near the phase transi- 
tion boundaries induced by either temperature or hole doping.0 The analysis of the exper- 
imental data showed the presence of mixed phases. For compositions intermediate near to 
the metal-insulator transition at low temperatures the ferromagnetic and antiferromagnetic 
phases coexist while the ferromagnetic and paramagnetic ones coexist in the vicinity of the 
paramagnetic to ferromagnetic transition temperature. On the other hands the existence of 
a homogeneous canted phase first proposed by de GenneS in the attempt of explaining the 
weak ferromagnetism of the phase with < x < 0.2 has not been confirmed. 

In these materials, the ferromagnetic phase is induced by the combined effect of the Cg 
electron hopping between nearest neighbor sites and Hund' s exchange with the localized 
t2g electron spins: this phase is associated with the properties of a metallic state. On the 
contrary the antiferromagnetism is induced by the exchange integral between the d shells of 
adjacent atoms and it is connected with the properties of an insulator state. 



In this paper we study, on the basis of a variational approach, the interplay of the 
electron-phonon interaction and of the double and super exchange magnetic effects as func- 
tion of the doping and the temperature in the range < x < 0.5. We will show that in 
the intermediate electron-phonon coupling regime the combined effect of the magnetic and 
electron-phonon interactions pushes the system towards a charge density instability. This 
explains the experimentally observed tendency of the manganites to form inhomogeneous 
magnetic structures near the phase boundaries due to temperature or Ca concentrations. 
The metal-insulator transition in these systems is a first order transitionlli accompanied by 
the coexistence of two phases: a high density one made by localized small polarons forming 
paramegnetic or antiferromagnetic domains depending on temperature and a low density 
one made by itinerant large polarons forming ferromagnetic domains. 

From a theoretical point of view the coexistence between hole-undoped and hole-rich 
phases has been discussed by using exact numerical approaches on small lattices treating 
the Jhan- Teller phonons classically^. However, little attention has been given to the quan- 
tum nature of the Jhan- Teller phonons, and, more important, to the possible formation 
of polarons characterizing the lattice deformations associated to the phase separation. On 
the other hand, after Millis et al.i, the role of polarons in the manganites has been clearly 
recognized. The main goal of the present letter is to show that the two scenario are fully 
compatible and that the formation of two types of polarons is responsible, together with the 
magnetic interactions, of the complexity of the manganites in the region < x < 0.5. 

The model. In Lai_xC a^MnO^, the electronically active orbitals are the Mn 3d orbitals 
and the mean number of the d electrons per Mn is 4 — x. Three electrons go into the core 
states and the remaining 1 — x electrons occupy the outer shell eg orbitals. The Jahn- Teller 
effect splits the e, double degeneracy^ and the single orbital approximation is reasonable 
for X < 0.5. The corresponding Hamiltonian is given by: 



H = -t 
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+ ^E4c.(a,-..-^' + at^--^'). (1) 

Here the first term describes the double exchange mechanism, in the limit where the 
intra-atomic exchange integral J is much greater than the transfer integral t, S'q"' represents 
the total spin of the subsystem consisting of two localized spins on sites i and j (nearest 
neighbors) and the conduction electron, 5* is the spin of the t2g core states, ojq denotes 
the frequency of the optical local phonon mode, the dimensionless parameter g indicates 
the strength of the electron-phonon interaction in the Holstein model,^ e represents the 
antiferromagnetic super-exchange coupling between two nearest neighbor t2g spins and is 
the number of lattice sites. The units are such that h = 1. 

The variational approach. 

We perform two successive canonical transformations. The first is a generalization of the 
Lang-Firsov unitary transformation^: 

where the phonon distribution function has to be determined variationally. 

This choice is related to the approximation first used by Tomonag a0 in the treatment 
of the coupling between mesons and nucleons and based on the assumption that all mesons 
are virtually emitted independently of the others. It recovers the well known Lang -FirsovS 
canonical transformation assuming narrow band small polaron, and, moreover, it provides, in 
the weak coupling regime, the generalization to the many-body problem of the unitary trans- 
formation utilized by Lee-Low-Pines0 in the continuum approximation within the Frohlich 
model.0 

The second canonical transformation is a Bogoliubov transformation. It introduces cor- 
relations between the emission of successive virtual phonons by the electron and it is re- 
sponsible of the phonon frequency renormalization^: 



5*1 = exp 



(2) 



5*2 = exp 



-a J2 - a-q^q) 



(3) 



where a indicates a variational parameter. 



The above introduced canonical transformations used to treat variationally the electron- 
phonon interaction are known to give a satisfactory description of the single polaron problem 
in all range of couplings being able to interpolate from weak (large polaron) to strong (small 
polaron) interactions. For this reason we adopt the same approach here for a system of 
interacting polarons. However, it is worth to note that, in the present case, the utilized 
approach is expected to be even better than in the single polaron problem. In fact, the choice 
of richer trial wave function made as a linear superposition of Bloch states translationally 
invariant that, in the single polaron problem, gives a highly accurate description of the 
polaron features in the intermediate electron-phonon coupling regime,0 does not provide any 
improvement in a many-polaron system since the asymptotic wave functions are orthogonal 
and the off-diagonal matrix elements of the Hamiltonian are zero in the thermodynamic 
limit. 

The transformed Hamiltonian is: 

/ 5"''-?' + 1 /2 \ 



^ c|q (a,- + al^-) {g - /,-) e^^'^^e^" + sinh 2a cosh 2a (a^^l^ + 



^—cflq 



+u^N (sinh 2af + ^ E (/| " 2^?/^) cos [q- [Ri - R, 



riirij + e'^Si- Sj 



(4) 



where 



XjXj = exp 



(5) 



By assuming a simple square of density of states, in the mean field approximation, the free 
energy per site becomes: 

F 



N 



±-ZS^ml + KbT log (l - e-'^'^«™^^4") + p (1 - p) 



N 



w 



+^0 



(sinh2a)^ - p'g' + KbT [p\ogp + (1 - p) log (1 - p) + Xnis - logz/5] (6) 



where the top and bottom sign, in the first term, hold, respectively, for the ferromagnetic and 
antiferromagnetic solutions for the localized spins. In the Eq. (0) Z indicates the number of 



nearest neighbors, p is the electron concentration, A is a dimensionless variational parameter 
that is proportional to the effective magnetic field, vs and ms represent, respectively, the 
partition function and the magnetization of the localized spins and W denotes the full band 
width: 

where a is the lattice parameter and Nq is the average phonon number with renormalized 
frequency uoq cosh 4a. The factor ( ( ^Is+i ) ) estimate of the reduction of the hopping 

due to the double exchange mechanism and it is evaluated in the ferromagnetic case by 
following Kubo and Ohatail and in the antiferromagnetic case by making use of the Clebsch- 
Gordan coefficients. Finally the expression of the phonon distribution function is obtained 
by minimizing the free energy per site: 

^ 1 + (2iVo + 1) [1 - cos (g.a)] ' 

The results. In the case of a single electron coupled to an optical local phonon mode 
we have shownill that a linear superposition of two Bloch wave functions describing the 
small and large polarons gives a very good estimate of the polaron ground state energy and 
provides an highly accurate description of the polaron features in the intermediate electron- 
phonon coupling regime and for electron and phonon energy scales not well separated. In 
particular, in the intermediate polaron the spectral weight is equally distributed between the 
ground and all the excited states and the polaron wave function contains, in equal weights, 
the small and large polaron contributions. It is naturally to ask what happens when there 
are many polarons since, now, we can expect a phase separation. To this end we first 
perform the minimization of the free energy, subject to the condition that the total number 
of particles is fixed, with respect to the two variational parameters a and A by making 
use of the self-consistent equation for the bandwidth W . Then, in order to investigate the 
existence of coexisting phases, we study the function G{p,p) = F{p) — pp, where F{p) 
is the free energy calculated from the above discussed minimization procedure and p is 



a Lagrangian multiplier to be adjusted to get the desired total number of particles. The 
function G{p,fi), at fixed value of the temperature, is reported in Fig.l for different values 
of the chemical potential /j, in the intermediate electron-phonon coupling regime Ap ~ 1 
(Ap represents the ratio between the small polaron binding energy and the energy gain of 
an itinerant electron on a rigid lattice). Prom Fig.l it is evident that there are values of 
the chemical potential such that the function G is stationary for three distinct electron 
densities. In particular, there is a value of /i such that G is the same at two different values 
of the electron density: pi and p2- They characterize the coexistence region. We stress that 
pi and p2 depend on the temperature and correspond to homogeneous phases of the system 
constituted by small and large polarons. Below a critical temperature the system segregates 
(see Fig. 2) in antiferromagnetic or paramagnetic and ferromagnetic domains (with density pi 
and P2) characterized by different values of the lattice distortions and by distinct behaviors 
of the resistivity with the temperature. The domains of large polarons are characterized by a 
weak electron mass renormalization and by coherent motion. On the other hand within the 
domains of small polarons the well-known band collapse takes place, the coherent motion 
is suppressed rapidly with increasing the temperature and the resistivity has a thermal 
activated behavior characteristic of the semiconducting phase. 

In the region of coexistence the fractions of volume, Vi/V and V2/V, filled with density 
Pi and P2, are determined by the two following conditions: Vi/V + V2/V — 1 and Vi/Vpi + 
V2/VP2 — P, where p = n/V and n represents the total number of charge carriers. It is 
crucial, in order to confirm our results, to evaluate the lattice displacement associated to 
the two kinds of coexisting domains. An estimate of the average deviation of the molecule 
on the site i from the equilibrium position when one electron is on site i is given by the 
function: 




(9) 



where M denotes the ionic mass. 



Applying the two canonical transformations 5*1 and 5*2 one finds: 
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D 



In particular, for narrow band small polarons is of order of g and D assumes the value 

By using the following values for the model parameters: t = 2uJo, g = 3.5, e = O.Olt, 
ujq = 50meV, for M the mass of the oxygen atom and by taking into account that the Jahn- 
Teller distortion of the Mn"^^ ion consists in an axial elongation of the two MnO bonds of 
the MuOq octahedra, at T = 2A0K we find: D/2 ~ 0.18 A for the paramagnetic domains 
of small polarons and D/2 ~ 0.12 A for the ferromagnetic domains of large polarons. These 
results are in agreement with the experimental data0 on LaQ -j^CaQ,2bMnO^ compound, that 
show evidence, above ITOi^', of small polarons characterized by MnO bond of order of 
~ 0.22 A and large polarons with MnO bond of order of ~ 0.1 A. 

In summary, we have studied the combined effect of the magnetic and electron-phonon 
interactions in Lai-xCa^MnO^ perovskites by means of two canonical transformations and 
the mean field theory. It has been shown that in the intermediate electron-phonon cou- 
pling regime these systems form inhomogeneous magnetic structures near the metal-insulator 
transition, segregating into antiferromagnetic or paramagnetic and ferromagnetic domains 
of small and large polarons characterized by different values of the lattice distortions as 
experimentally observed. 



FIGURE CAPTIONS 

Fl The function G{p,fi) = F{p) — pp is plotted as a function of the doping for different 
values of the chemical potential at fixed temperature: KbT = OAuq (for uq = 50meV 
KbT = O.Scuq corresponds to a temperature T of order of SOOi^). 

F2 The phase diagram corresponding to: t = 2ujo, g = 3.5, e = O.Olt. The area between the 
dotted lines indicates the region where the localized and delocalized phases coexist. 



8 



^G.H. Jonker and J.H. Santen, Physica 16, 337 (1950); J.H. Santen and G.H. Jonker, Physica 
16, 599 (1950). 

2 P. SchifFer, A.P. Ramirez, W. Bao and S.W. Cheong, Phys. Rev. Lett. 75, 3336 (1995). 

^C. Zener, Phys. Rev. 81, 440 (1951); C. Zener, Phys. Rev. 82, 403 (1951). 

^A.J. MiUis, P.B. Littlewood and B.I. Shraiman, Phys. Rev. lett. 74, 5144 (1995); A.J. Milhs, 
B.I. Shraiman and R. Mueller, Phys. Rev. Lett. 77, 175 (1996). 

5 P.W. Anderson and H. Hasegawa, Phys. Rev. 100, 675 (1955); P.G. de Gennes, Phys. Rev. 118, 
141 (1960). 

6 Guo-Meng-Zhao, K. Coder, H. Keller and K.A. Muller, Nature (London) 381, 676 (1996). 

^A. Lanzara, N.L. Saini, M. Brunelli, F. Natali, A. Bianconi, P.G. Radaelli and S.W. Cheong, 
Phys. Rev. Lett. 81, 878 (1998). 

^K.H. Kim, J.Y. Gu, H.S. Choi, G.W. Park and T.W. Noh, Phys. Rev. Lett. 77, 1877 (1996). 

9 A.J. Millis, Nature (London) 392, 147 (1998). 

°G. Allodi, R. De Renzi, G. Guidi, F. Licci, and M.W. Pieper, Phys. Rev. B 56, 6036 (1997); 
G. Allodi, R. De Renzi, and G. Guidi, Phys. Rev. B 57, 1024 (1998); J. Dho, I. Kim, S. Lee, 
K.H. Kim, H.J. Lee, J.H. Jung, and T.W. Noh, Phys. Rev. B 59, 492 (1999); S.J.L. Billinge, 
Th. Proffen, V. Petkov, J.L. Sarrao, and S. Kycia, Cond-Mat 9907329 (1999). 

^ P.G. de Gennes, Phys. Rev. 118, 141 (1960). 

^ J. B. Goodenough and J. S. Zhou, Nature 386, 229 (1997). 

3 A. Moreo, S. Yunoki and E. Dagotto, Science 283, 2034 (1999) 

^ J. Zhang, A.R. Bishop, and H. Roder, Phys. Rev. B 53, 8840 (1996). 

^T. Holstein, Ann. Phys. 8, 325 (1959); T. Holstein, Ann. Phys. 8, 343 (1959). 



9 



I. J. Lang and Yu. A. Firsov, Soviet Physics JETP 16, 1301 (1963); Yu. A. Firsov, Polarons 
(Moskow, Nauka, 1975). 

S. Tomonaga, Prog. Theor. Phys. 2, 6 (1947). 

T.D. Lee, F. Low, and D. Pines, Phys. Rev. 90, 297 (1953). 

H. Frohhch et al., Philos. Mag. 41, 221 (1950); H. Frohhch, in Polarons and Excitons, edited by 
C.G. Kuper and G.A. Whitfield (Ohver and Boyd, Edinburg, 1963), p. 1. 

H. Zheng, Phys. Rev. B 36, 8736 (1987); H. Zheng, Sohd State Communication 65, 731 (1988). 

V. Cataudclla, G. Dc Filippis, and G. ladonisi, Phys. Rev. B 60, 15163 (1999); V. Cataudella, 
G. De Filippis and G. ladonisi, Phys. Rev. B 62, 1496 (2000). 

K. Kubo and N. Ohata, J. Phys. Soc. Jpn 33, 21 (1972). 



10 



T 



IT) 



I 

II 




X 



d 



d 



I 



I I I I I 




CD 



00 



I 



I 



CNj 

T 



Csl 



to 



00 
C\l 




X 



CD(N^t£>0Or\(N^CD 

' T T T T ' T T T 




T 




X 



in 

-I d 



in 



X 



-1.46 
-1.48 


-1.5 
-1.52 
-1.54 
-1.56 






in 




T 




T 




II 






iTi 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



X 



d 



d 



^ Csl 

I p 



o 



o 



00 ^ 

o 



incDini^LnooincDin^ 
mmcDroi^moorocn 

T T T T 

I I I I I 



00 

d 



CL 



6 



CD 

6 



in 
d 



d 



+ 

CL 



d 




d 



o 



